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ABSTRACT 

The recent observation of the pulsar PSR J1614-2230 with a mass of 1.97 ± 0.04 M0 gives a strong 
constraint on the quark and nuclear matter equations of state (EoS). We explore the parameter ranges 
for a parameterized EoS for quark stars. We find that strange stars, made of absolutely stable strange 
quark matter, comply with the new constraint only if effects from the strong coupling constant and 
color-superconductivity are taken into account. Hybrid stars, compact stars with a quark matter core 
and an hadronic outer layer, can be as massive as 2 M©, but only for a significantly limited range of 
parameters. We demonstrate that the appearance of quark matter in massive stars depends crucially 
on the stiffness of the nuclear matter EoS. We show that the masses of hybrid stars stay below the 
ones of hadronic and pure quark stars, due to the softening of the EoS at the quark-hadron phase 
transition. 

Subject headings: stars: neutron — equation of state 



1. INTRODUCTION 

The densities in the interior of neutron stars exceed the 
ground state density of atomic nuclei, no ~ 0.16 fm~'^, 
by far. This naturally raises the idea, that compact stars 
might contain a deconfined and chirall y restored quark 
phase. Recently, iDemorest et al.l (|2010( ) found a new ro- 
bust mass limit for compact stars by determining the 
mass of the millisecond pulsar PSR J1614-2230 to be 
M = 1.97 ± 0.04 Mq. This value, together with the 
mass of pulsar J1 903-H0327 of M = 1.667 ± 0.021 Mq 
(jFreire et al.ll2010t ) is m uch larger than the Hulse-Taylo r 
limit of M ~ 1.44 M© (jThorsett fc Chakrabart"vlll999f ). 
which for a long time has been the highest precisely mea- 
sured pulsar mass. In this letter we want to explore the 
implications of this new measurement on the possible 
presence of quark matter in compact stars. Moreover, 
our aim is to map out the parameter range for the widely 
used quark bag model with respect to its ability to repro- 
duce high mass compact stars such as PSR J1614-2230. 
There are two classes of compact stars which contain 
quark matter. The first class are so-called hybrid stars, 
with quarks only in their interior either in form of a 
pure quark matter core or a quark-hadron mixed phase. 
The size of the core depends hereby on the critical den- 
sity for the quark-hadron phase transition ricrit under 
neutron star conditions. The second class of so-called 
(strange) quark stars is realized for the special sce nario 
of abs o lutely stable str a nge qua r k mat ter (see e.g. Iltohl 
([T970I) : IBodmeii (119711) : IWittenI (IT98l 'l. It is based on 
the idea that the presence of strange quarks can lower 
the energy per baryon of the mixture of up, down, and 
strange quarks in weak equilibrium below the one of 
^^Fe (~ 930 MeV). As a consequence, this strange quark 
matter forms the true ground state of nuclear matter 
and occupies the en tire compact star (jAlcock et al.lll986l : 
IHaensel et aIIIT986l ). 

The mass measurement for PSR J1614-2230 sets for 
the first time very strong limits for the parameters 
of any zero temperature equation of state (EoS), and 



thereby also for the one of quark matter. Usually, 
the appearance of strangeness in quark and hadronic 
matter provides an additional degree of freedom and 
thereby softens the nuclear EoS, that is, decreases 
the pressure for a given energy density. As a re- 
sult, quark and hybrid stars cannot reach high masses. 
However, many studies found that effects from the 
strong interaction, such as one-gluon exchange or color- 
superconductivity can stiffen the quark matter EoS 
and i ncrease the maximum mass of quark and hybrid 
stars fRiister fc Rischke"200 4l: iHorvath fc Lugoneij |2004l: 
[Affor d ct al. 2007; Fis cher et al.l 120101: KurkclaetaLl 
IMOallU) . lOzel et all (|2010l) and iLattimer fc PrakashI 
(|2010D gave first studies on the implications of the new 
mass limits from PSR J1614-2230 for quark and hybrid 
stars in the quark bag model. However, as we will show 
below, a systematic analysis of the whole allowed param- 
eter range is still missing. 

ILattimer fc PrakashI ()2010[ ) include strange quark mat- 
ter in form of a bag model EoS for quark stars as well 
as hybrid stars. The authors do not study strong ef- 
fects from color-superconductivity and impose the ad- 
ditional constraint of 
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Lattimer fc Pra kash 2010;), but not necessary condition 
( Wittenlll984D . Moreover, they exclude a priori the ex- 
istence of a quark-hadron mixed phase and come to the 
conclusion, that the existence of a 2.5 M© star would ex- 
clude the quark-hadron phase transition in compact star 
interiors. As we will show in the next sections our results 
cannot confirm this statement, furthermore we find that 
a quark-hadron mixed phase in fact plays a major role 
in supporting high mass hybrid stars. 
In a different analysis by lOzel et aD ()2010( ) , the authors 
studied the implications of the new measurement on hy- 
brid stars with a parameterized quark bag model includ- 
ing effects from color-superconductivity and QCD cor- 
rections. They find that both effects are required to sup- 
port the mass of PSR J1614-2230. However, lOzel et al.l 
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(|2010f ) adjust the bag constant to obtain a fixed density 
of ricrit = l-5no for the phase transition to quark matter 
from the relatively soft APR nucleonic EoS (|Akmal et alj 
Il998jl. We find that the stiffness of the hadronic EoS is 
important for large hybrid star masses and also that the 
maximum mass of hybrid star configurations experiences 
a minimum at around ricrit = 0.1 — 0.2fm~'^ - values close 
to the critical density which [Ozel et alJ (|2QlClD choose for 
their calculations. 

Therefore, the aim of this paper is to fully and systemat- 
ically exploit the constraints on the quark bag EoS pro - 
vided by the new mass limit of lDemorest et all (j2010[ ). 
Quark matter is described by a bag model EoS with 
first order corrections from the strong interaction cou- 
pling constant and effects from finite strange quark mass 
and color-superconductivity. For the hybrid star calcu- 
lations, we use the two diff erent relativi stic mean-field 
(RM F) parameter sets TMl (jSugahara fc Toki.1994. ) and 
NL3 (jLalazissis et al.|[l997D . to explore the influence of 
the hadronic part of the EoS. In our calculations we do 
not include hype rons which can alter the qu ark hadron 
phase transition (jBhattacharvva et al.ll20l"o() . However, 
their exact role is cu rrently an open question (see e.g. 
lYasutake et al.l ()20ICl[ )). Therefore, in this work, we will 
focus on non-strange hadronic matter. We consider the 
two possible extreme cases for the phase coexistence be- 
tween quark and hadronic matter: the Maxwell tran- 
sition, corresponding to a very large surface tension of 
quark matter (Heisclberg et al.; |1993| ). and the Gibbs 
construction ( Glendenning. ,19921 ) which completely ne- 
glects Coulomb and surface energies. 
In the following we will describe our results and compare 
them with the aforementioned studies. Sections 12.11 and 
12.21 are devoted to quark stars with unpaired and color- 
superconducting quark matter in the color-flavor-locked 
(CFL) phase, respectively. Hybrid stars are discussed in 
section [3l 

2. QUARK STARS 

2.1. Unpaired quark matter 

For the strange quark matter, we take the modified 
bag model: 

o 4 

^QM= ^^ + £l(l-«4)+Seif, (1) 

where fi^ are the Grand potentials for the up, down, 
and strange quarks and electrons describing these as 
non-interacting fermions. We choose the stra nge quark 
mass to be = 100 MeV (jAmsler et al.ll2008D while the 
masses of the up and down quarks and electrons are set to 
zero. In the s ense o f the generic quark matter EoS from 
lAlford et all (gOOS^, we have added the 04 term with 
the baryon chemical potential /i of the quarks in order 
to account for corrections from strong interaction. The 
usual approach in quark bag models is to unite all non- 
perturbative effects of the strong interactions into a bag 
constant B. The EoS can then be extended by including 
first order corrections in the strong coupling constant (see 
e.g. Fragaetal. 2001). The quark bag model in equation 
((1} is motivated by this approach. However, since quark 
star matter is not in the perturbative regime, we consider 
04 and the bag constant as effective parameters, denot- 
ing the latter by Bf,ff , and explore their whole parameter 
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Fig. 1. — Maximum masses of unpaired strange quark stars as a 
function of B^g and 04. The green shaded area marks the allowed 
parameter region according to the constraints of the existence of 
nuclei (2-flavor line), absolute stability of strange quark matter (3- 
flavor line), stability of fast rotating stars (Kepler line), and the 
mass of PSR J1614-2230 including its la error. 

range. Therefore, we vary 04 from 04 = 1, which corre- 
sponds to no QCD corrections, to small values when the 
corrections are strong. Equation ([T]) enables us to com- 
pute the pressure, energy density, and baryon number 
density assuming charge neutrality and /3— equilibrium. 
By solving the Tolman-Oppenheimer-Volkoff equations 
we obtain the maximum quark star masses. 
Following iFarhi fc Jaffi (Il98l . we require nonstrange 
quark matter in bulk to have a binding energy per baryon 
higher than the one of the most stable atomic nucleus, 
^®Fe, which is 930 MeV, plus a 4 MeV correction coming 
from surface effects. By imposing that E/A> 934 MeV 
for 2— flavor quark matter at ground state, we ensure 
that atomic nuclei do not dissolve into their constituent 
quarks. Thereby we obtain an upper limit on the max- 
imum mass of strange quark stars denoted as "2-flavor 
line" in Fig. [TJ Another constraint is given by the im- 
plementation of the strange matter hypothesis (jBodmed 
IT97I; Wittcn 1984) as described in the introduction, with 
E/A < 930 MeV f or str ange quark matter at ground 
state (Farhi fc Jaffdl[T98^ . This condition results in the 
"3- flavor line" in Fig. [1] and gives a lower limit on the 
maximum masses. Fig. [1] also shows lines of constant 
maximum mass. The three dotted lines enclosing the 
red shaded area re present the mas s of PSR J1614-2230 
with its la error (jPemorest et al.l 12010 ). The 2-flavor 
and the 3-flavor lines cross on the left outside the plot 
range at 04 = 0.247, Slj^ = 102.24 MeV which corre- 
spond to a maximum mass star with M= 3.36 Mq and 
a radius of 19 km. The Kepler line at low B^jj repre- 
sents a limit for quark stars which can rotate with a 
Kepl erian frequency of at least 716 Hz (jHessels et al.l 
I2006f) . Therefore, the green shaded area is the allowed 
quark star parameter region with a maximum mass of 

2.54 Mq at 04 « 0.53 and bI^j^^ « 123.7 MeV. However, 
the Ke pler line is obtained from a parametrization of 
ILattimer fc Prakas h (2 007.) and gives a rough estimate 
when applied to strange stars. For a more reliable Ke- 
pler limit, the presented quark EoSs should be applied 
in general relativistic c alculations of rotatin g q uark stars 
similar to the studies of^ Haensel et al.l ()2009[ ) or lLo fc LinI 
(|2011i) . From Fig[Tl we see that for 04 = 1 the 2-flavor 



line requires M^ax ^ 1-92 Mq which is ruled out by the 
new mass limit, at least within its la error. Thus we find 
that 04 < 1, i.e. QCD corrections must be included to 
ensure the compatibility of the model with observational 
data. 

2.2. Color- superconducting strange matter stars 

At large densities, such as in compact star interiors, 
up, down, and strange quarks are assumed to undergo 
pairing and form the so-called C FL phase. We adopt 
the EoS from Aff ord et al.l ()2001[ ) which introduces the 
pairing energy A as a new free parameter: 

6 r 
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where v = 2p, — ^ ^j? — m^/S. We added again the 
04— term to account for QCD corrections. The results 
for 04 — 1.0 and = 100 MeV are shown in Fig. [5] 
where we have imposed the same constraints as in Fig. [TJ 
The " 3- flavor" and " 2-flavor" lines give again a lower and 
upper limit on the maximum mass. The green shaded 
area is the maximally allowed parameter region. Its up- 
per edge is solely given by the constraint from the mass 
measurement of PSR J1614-2230. Note, that for small 
values of A quark stars are not allowed. For large gaps 
starting at A > 20 MeV the allowed area opens up and 
one can obtain high maximum masses. The largest mass 
allowed within the plot range is 2.34 M© at 04 = 1.0, 

Bl'g = 145 MeV and A = 100 MeV. If we assume the 
same constraint from the Keplerian frequency as before, 
we find, that it has no influence on our results within the 
plot range of Fig. [2] More exotic parameter combina- 
tions, as e.g. 04 = 0.66, B\'g = 130.5 MeV, A = 50 MeV 

and 04 = 0.75, B\'g = 134.9 MeV, A = 100 McV with 
maximum masses of 2.5 M© and 2.8 M© respectively are 
also allowed. 

Taking into account QCD corrections, i.e lowering 04, 
gives maximum mass lines shifted to only slightly lower 
values of B^g- Together with the lowered 2-flavor and 
3-flavor lines this means that at some point the 3-flavor 
constraint will become important. Still we can obtain 
quite large maximum masses at a sufhciently high gap 
value. Varying the strange quark mass basically results 
in shifting the maximum mass lines along the gap-axis as 
the crucial contribution to th e EoS comes from a term 
ml - 4A2 (I Afford et al. l [2001 . Thus, for fixed 04 and 
a higher strange mass, the allowed area opens up at a 
larger value of A. 

3. HYBRID STARS 

For the hybrid star calculations we use again the bag 
model EoS of equation ^ and the same quark masses 
as in the previous sections. As hadronic EoSs we choose 
the TMl and NL3 RMF parameter sets, with maximum 
neutron star m asses of 2.2 M© a nd 2.78 M©, r espectively 
(jSugahara fc T oki 1994; Lalazis sis et al.lll997t ). The crit- 
ical density Ucrit for the quark-hadron phase transition 
is mainly given by the effective bag constant Bf.g while 
corrections from the strong interaction 04 affect the stiff- 
ness of the quark EoS. Combinations of small values for 




Fig. 2. — Same as in Fig. [T] but this time for color- 
superconducting strange quark matter with the pairing gap A and 
fixed a4 = 1.0. 

04 and Bf.ff lead therefore to stiff hy brid EoSs with a low 
critical density (jFischer et al.|[20lol ). As a consequence, 
the corresponding hybrid stars can be very massive and 
have a large quark matter core. Similar to lOzel et "aP 
(|2010D we omit all parameter combinations which lead 
to several transitions between quark and hadronic mat- 
ter, which happens for small B^^g and 04. 
The phase transitions are modell ed by the Gi bbs and 
Maxwell constructions (see e.g. IGlendennind ([1992)) 
corresponding to a small (Gibbs) and to a large value 
of the surface tens ion (Maxwell) (see discussion in 
iPage fc Reddvl (|2006l) ). The former leads to the presence 
of an extended quark-hadron mixed phase region in the 
compact star interior where the pressure rises smoothly 
with density. With the Maxwell approach, matter experi- 
ences a direct transition from hadronic to quark matter in 
cold compact stars, accompanied by a density jump from 
lower (hadronic phase) to higher (quark phase) densities. 
Applying both models for the phase transition we calcu- 
late the hybrid star maximum masses and plot them as 
a function of the bag constant B^g for fixed values of 
04. Fig. [3] shows the hybrid star maximum mass curves 
for the TMl EoS for the hadronic phase, while in Fig. |4] 
we show the results for the NL3 model. The lines which 
extend from low values to high values of Sejj correspond 
to calculations using the Gibbs phase transition. Stars 
on the solid lines have a pure quark matter core while the 
dashed lines represent stars where only a mixed phase is 
present. The maximum masses for the Maxwell transi- 
tion are represented by the grey shaded area. Due to the 
absence of a mixed phase, hybrid stars in the Maxwell 
approach can only contain a pure quark matter core. A 
too large density jump from hadronic to quark matter in 
their interior leads to a gravitational instability against 
radial oscillations. As a consequence we find from Figs. 
[3] and m that the parameter range for hybrid stars in the 
Maxwell approach is significantly reduced in comparison 
to the "Gibbs hybrid stars". Furthermore, in most of the 
cases we find that for the same combinations of B^g and 
04 the Maxwell phase transition leads to lower hybrid 
star maximum masses than the ones with a Gibbs con- 
struction. Only for low values of 04 when quark matter 
becomes very stiff, the maximum masses of the " Maxwell 
hybrid stars" can become significantly larger than their 
Gibbs counterparts. 
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Fig. 3. — Maximum masses and critical densities of hybrid stars 
calculated with the TMl RMF EoS as a function of B^ff and 04. 
The quark-hadron phase transitions are modelled by the Maxwell 
and Gibbs approach. For the latter, solid lines indicate pure quark 
matter cores and dashed lines indicate mixed phase cores. 




.0 .0 Jg 

34=1 .0 



» 



PSR J1614-2230 



Maxwell hybrid stars 

Maxwell, quark core 

Gibbs, quark core 

Gibbs, mixed core 

~ 0.1 1/fm, ♦ 
n,~ 0.2 1/fm, O 
n^- 0.3 1/fm, • 
n.~0.4 1/fm^ O 
nr~0.5 1/fm'' A 



200 220 
Be,,"" [MeV] 



240 



260 



280 



Fig. 4. — Same as in Fig. |3]with hadronic matter described by 
the NL3 RMF EoS. 



For stable hybrid stars with the Maxwell and Gibbs tran- 
sition, the pure quark and the mixed phase extend over 
almost the whole star when B^ff, i.e. the critical den- 
sity ricrit, is low. However, as can be seen from Figs. 
[3] and m a pure quark matter core in hybrid stars oc- 
curs only for small values of Bf^g. For the TMl hadronic 
EoS, sufficiently massive Gibbs hybrid stars contain only 
a mixed phase in the core. On the other hand, for the 
NL3 hadronic EoS with low values of B^g and Ucrit , large 
cores of pure quark matter exist down to 04 ~ 0.5. This 
is a consequence of the much stiffer NL3 hadronic EoS in 
comparison to the TMl parameter set. Nevertheless, hy- 
brid stars with a pure quark core and a mass > 1.93 Mq 
are only obtained for 04 < 0.6. 

A common feature which is seen for hybrid stars in the 
Gibbs approach is that for a fixed value of 04 the max- 
imum masses decrease with lower critical densities and 
experience a minimum around ricrit ~ 0.2 fm^"^. For 
smaller ricrit, the quark EoS starts to dominate and the 
maximum masses increase again as they approach the 
limit of absolutely stable strange quark matter. We plot 
the mass of PSR J1614-2230 with its la error as a gray 
band as weU as lines at M= 1.44 M0 and M= 1.67 M© 



to indicate the masses of the Hulse- Taylor pulsar and of 
J1903-h0327 respectively (jThorsett fc Chakrabartvl 19991: 
lFreirel[2009l ). Especially in Fig. [3] it can be seen that 
the new mass limit significantly tightens the constraints 
on the model parameters, as the whole area below the 
gray band is now excluded while formerly this was only 
the case below the 1.67 Mq line. In this work we did 
not consider effects of a finite surface tension on the 
quark-hadron mixed phase (jHeiselberg et al.lll993() which 
has been found to be an intermediate of the Gibbs 
and the Maxwell constructions ( Maruva ma et al.l 120071 : 
lEndo et al.l 1200(1 ). However, the B^ff — parameter 
space for stable hybrid stars with finite surface ten- 
sion can be expected to be between the Gibbs and the 
Maxwell constructions. Therefore our result, that hybrid 
stars can be massive, remains valid. 

4. CONCLUSION 

We present for the first time a comprehensive and 
systematic study on the constraints of the new com- 
pact star mass limit from the millisecond pulsar PSR 
J1614-2230 on the properties of quark and hybrid stars 
modelled within an extended quark bag model. The 
parameters of the bag model are an effective bag con- 
stant, corrections from the strong interaction coupling 
constant and color-superconductivity. We find that 
the new mass limit does not rule out the possibility of 
having quark matter in compact stars but provides tight 
bounds on its properties. High compact star masses 
where quark matter is the dominant component, require 
strong QCD corrections and/or a large contribution 
from color-superconductivity. In this case strange stars 
can reach masses far beyond 2 Mq. For hybrid stars 
with a sizeable quark matter phase, our investigation 
shows that pure quark matter cores are obtained only 
for a small parameter range when the hadronic EoS is 
stiff and the critical density for the quark-hadron phase 
transition is around saturation density. Our results 
agree with lOzel etHI (pOlOl ) concerning the importance 
of effective QCD corrections to reach high compact 
star masses. C ontrarily to one of the statements of 
iLattimer fc Prak ash (2 010l ). we find that pairing helps 
to increase the maximum mass and that corrections 
from the strong interaction have a significant effect. For 
hybrid stars we demonstrate that the allowed parameter 
region hinges crucially on the stiffness of the hadronic 
EoS and can therefore be much larger than in the case of 
lOzel et all (|2010D . Recently, Ivan Kerkwijk et al.l (|2011l ) 
reported the possible existence of a 2.4 Mq star. If such 
a measurement is confirmed, even larger corrections 
from the strong coupling constant and larger values of 
the CFL gap will be required for strange stars. Hybrid 
stars could exist only for a stiff hadronic EoS and would 
contain only a core with a quark-hadron mixed phase in 
our approach. 

An investigation similar to the one we presented here 
would also be desirable for other effective models of QCD 
such as the Nambu -Jona-Lasini o (NJL) and the PNJL 
models (iKlahn et al. 2007; Paeliara fc SchafFner-BielichI 
I200I iBlaschke et al.lHoiS llppolito et al.ll2008D . as weh 
as Schwinger-Dyso n approaches, as recently proposed 
bv lLi et al.l (I2OIII) . Although mass measurements are 
very useful for constraining the nuclear matter EoS, 
additional observational information is required to probe 
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the existence of quark matter in compact stars. Cooling, 
r-modes cal culations, and gravi tational wave signals 
of mergers (|Bauswein et alj I2OO90 are promising tools. 
Also heavy ions collisions experiments provide crucial 
information on the nuclear matter EoS: an extended 
analysis on the quark models paramaters which includes 
both astrophysical constraints, as the mass of PSR 
J1614-2230, and terrestrial laboratories constraints 
would be extemely interesting. 
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